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Abstract 


This computer subroutine computes the visual, blue and ultraviolet photometric 
magnitudes of the planets, their natural satellites, and the sun at varying observa- 
tion distances and phase angles. Currently available observational magnitude and 
phase function data are stored in the program and used in the computation. 
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UBV: Subroutine to Compute Photometric Magnitudes of the 

Planets and Their Satellites 


I. Introduction 

Design and performance evaluation of optical sensors 
and instruments for planetary missions requires the 
knowledge of planet and moon brightness. A computer 
subroutine has been written to calculate the Johnson and 
Morgan (Ref. 1) visual (V), blue (B), and ultraviolet (U) 
photometric magnitudes of the planets, all their known 
natural satellites, and the sun. 

II. Computation of Visual Magnitude 

The visual magnitude of a planet or satellite is given 
by: 

V = V(l, 0) + 5 log r^To - 2.5 log (<j>) (1) 

where V(l, 0) is the visual magnitude at one astronomical 
unit (1 AU = 149, 597, 893 km) from the sun and ob- 
server and 0 deg phase angle, r,, is the sun-body distance 


in AU, To is the observer-body distance in AU, and is 
the visual phase function dependent on the phase angle 4> 
where 0 < <j> < 180 deg. The phase angle is the sun- 
body-observer angle. All logs are to the base 10. Figure 1 
illustrates the observational geometry. 

For the sun the visual magnitude is given by: 

V = V,, + 5 log To (2) 

where V(l, 0) is replaced by V* the visual magnitude of 
the sun at one AU and To is the sun-observer distance 
in AU. 

The above equations are valid only for large planet- 
observer distances, typically greater than twenty planet 
radii. For distances closer than this, the phase angle 
varies over the surface of the planet and an integration 
technique is required. 
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PLANET OR 
SATELLITE 


HI. Stored Data 

Data stored in the program are listed in Tables 1 and 2. 
The data were obtained from various references and 
these are indicated by the numbers in parentheses in 
Table 1 and in the column headings of Table 2. 

Table 1 lists the object, its body number as assigned 
in the program, V(l, 0), the B-V and U-B color indexes, 
and the visual phase function assumed as identified 
by Table 2. The values for V(l, 0) are average values 
since actual magnitude can vary with the variation in 
the sun’s magnitude or with changing features due to 
planet or moon rotation, seasonal variations, atmospheric 
conditions, or other surface changes. This variation can 
be as large as two magnitudes as is the case with lapetus 
(Ref. 2). The value for Saturn is for “no-rings” as the 
actual magnitude varies with the changing aspect of 
the rings. An empirical formula (Ref. 2) for the change 
in magnitude due to the rings is 

AV = -2.60 sin B + 1.25 sin^ B {0° < B < 27”) (3) 

where B is the saturnicentric latitude of the observer. 

Additional variations in the values for V(l, 0) are pos- 
sible due to observational uncertainties. Since not all the 
planets can be observed at exactly 0 deg phase angle, 
some uncertainty is possible when extrapolating to 0 deg. 
Some of the dimmer moons have only been observed 
photographically, and errors may be introduced when 
converting to photometric magnitudes. 

Note that constant values of the color indexes B-V and 
U-B are stored in the program and that no computation 


is done. Although the indexes do vary with phase, there 
are only limited data available, and thus it is not mod- 
eled in the program at this time. Where no color infor- 
mation was available, a reasonable value has been stored 
in the program. These assumed values are noted by 
parentheses in Table 1. Uneertainties in the color indexes 
are due to the same factors that affect V(l, 0). 

Table 2 indicates how the visual magnitude varies 
with phase angle. The quantity —2.5 log (^) is also 
given by Am, and it is this quantity that is stored in the 
program. The phase function data have been determined 
empirically from observational data or from theoretical 
considerations. For planets and moons outside the earth’s 
orbit, observations can only be made over a limited 
range in phase angle. Mars can only be observed from 
the earth up to about 48 deg phase angle, and the value 
of Am must be extrapolated as indicated by the paren- 
theses in Table 2. Since the phase angle of the outer 
planets and their moons never exceed more than a few 
degrees as seen from the earth, a theoretical phase func- 
tion must be assumed. The same phase function has been 
assumed (Ref. 2) for the larger planets with optically 
thick atmospheres; Jupiter, Saturn, Uranus and Neptune. 
Since the size of Pluto is closer to the size of Mars, the 
Mars phase function has been assumed for Pluto. Since 
the phase function of the earth’s moon is the only one 
that has been accurately observed over all phase angles, 
it has been assumed as the phase function for all the 
moons at this time. 

Extrapolated values of Am at 180 deg phase, in Table 2, 
are given to allow interpolation in the program. The 
actual magnitude of a planet or moon at 180 deg phase 
would be a function of the observer-body distance, the 
sun-body distance, the planet or moon radius and the 
atmospheric properties of the planet or moon. These 
factors would determine if a halo were visible or not. 
Thus, care should be used when interpreting the pro- 
gram output when the phase angle is within a few de- 
grees of 180 deg. 

As a matter of interest Am is plotted for all the planets 
and moons in Fig. 2. 


IV. Program Usage and Logic 

The subroutine is coded in Fortran V for use on 
the JPL Scientific Computing Facility Univac 1108 
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Table 1. Planet and satellite body numbers and 
photometry data 


Color indexes 



Body 

number 


III. Titania 

IV. Oberon 

V. Miranda 

8. Neptune 

I. Triton 

II. Nereid 

9. Pluto 



Visual 

Color indexes 

V(1,0) 

B-V 

U-B 

+ 1.30(2) 

0.62 

0.25 (2) 

+ 1.49(2) 

0.65 

0.24 (2) 

+ 3.8 (2) 

(0.8) 

(0.4) 

- 6.87(2) 

0.41 

0.21 (2) 

- 1.16(2) 

0.77 

0.40 (2) 

+ 4.0 (2) 

(0.8) 

(0.4) 

- 1 .01 (2) 

0.80 

0.27 (2) 


Table 2. Values for Am = —2.5 log ((/>) 


Phase functions^ 

1 2 3 4 5 6 

(Ref. 3) (Ref. 3) (Ref. 3) (Ref. 5) (Ref. 2) (Ref. 3) 


^The numbers in parentheses indicate 
which the data were obtained. 


the various references 




^The column-head numbers refer to the planets and their moons, as 
follows: 

1 Mercury 4 Mars and Pluto 

2 Venus 5 Jupiter, Saturn, Uranus, and Neptune 

3 Earth 6 Moons 
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Computer. The subroutine is called by: CALL UBV 
(BN, PHA, RS, RO, V, BV, UB) where the input is: 

BN = body number as given in Table 1 

PHA = ^ phase angle of body in degrees, where 
0 < <j> < 180 deg 

RS = r, sun-body distance in AU 

RO = r„ observer-body distance in AU 

and the output is: 

V = the visual magnitude as given by equation 1 
or 2 

BV = B-V color index 
UB = U-B color index. 


The subroutine first checks the body number to see 
if it is acceptable. When an incorrect body number is 
input, the message “ILLEGAL BODY NUMBER” is 
printed, and the following values are output: 

V = 99.0 
BV = 0.0 
UB = 0.0 

When a correct body number is input, the program 
extracts the correct values of V(l, 0), B-V and U-B as 
shown in Table 1. V is computed using equations 1 or 2. 
The function Am is determined by linear interpolation 
between the values given in Table 2. 



PHASE ANGLE, deg 

Fig. 2. Variation of magnitude with phase angle for the 
planets and their moons 
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5UBR0UT I NE UBV[ BN , PHH , R5 , RD , V , BV , UB ! 


C VERSION OF 11/12/70 


DIMENSION NBr42) ,Rf 3,42) ,PHIf 6,19) ,NP(42) 


C R[1 ,J) *5 RRE Vfl ,0) ’5 


DRTH CRH ,J) .Jrl ,421/-2B.74,-.42,-4.41 , -3 .37 , . 21 , -1 . 73 , 12 ■ 1 , 1 3 . 3 , 
-9.25.-1 .90,-1 .53,-2.16,-! .20,6.30,7.00,9,30,12.1 ,11 .6,11 ,9,11 .4, 
12.1 ,-3.33,2.60,2.22, ,72, ,B9, .21 ,-l .16,4.61 ,1 .43,4 .45,4.45,-7.19, 
1.70,2,60,1.30,1.49,3.80.-6.37,-1.16,4.00,-1.01/ 

I ^ 1 

IC ^RC2,J)'5 RRE B-V’3 



DRTR [ RC 2, J) , Jr 1 ,42)/ .63, ,97, .82, -20. .92,1 .33, .60, ,60, .83,1 .17, 
.37, .33, .86, -SO. .80, .30, -30, .80, ,80, .SO, .30,1 .04, .60, .62, .73, .71 , 
.76,1 .30, .69, .71 , .80, ,80, .58, .80, ,80, .62, -65, .80, .41 , . 11 , .80, .80/ 


Rf3, J) *5 RRE U-B’5 


ORTH 

[RC3, J) 

, J=1 

,42)/. 14, 

, .40 

, .50, 

.00, 

, .46 

, .63, 

, .30, 

.30, 

.48 

,1,30, .52, 

.50, 

.55, .40, 

.40, 

.40, .40, 

.40, 

.40, . 

40, 

.40. 

.58, 

.30, . 

30, . 

34 , 

.30, .35, 

.75, 

.42, .28, 

.40, 

.40, .28, 

.40, 

.40, . 

25, 

.24, 

. 40 I 

.21 , . 

40 I . 

40, 

.27/ 


C; NP'S RRE MODEL NUMBERS 


DRTfl CNPCJi,Jrl, 421/0,1,2, 3, 6, 4, 6, 6, 5, 6,6,6, 6,6, 6,6, 6, 6,6, 6.6, 5,6, 

6,6,6,6,6,6,6,6,6,5,616,6,6,6,5,6,6,4/ 


NB'S CORRESPOND TO PLANET NUMBERS 



DRTR NB/0, -100, -200, -300, -301 ,-400,-401 ,-402,-500,-501 ,-602,-503, 
-504,-505,-506,-507,-508,-509,-510,-511 ,-512,-600,-601 ,-602,-603, 
-604 , -605 , -606 , -607 , -608 , -609 ,-610, -700 , -701 , -702 , -703 , -704 , -705 , 
-800,-801 ,-802,-900/ 


C PHI 'S RRE -2.5*L0G10CPH] ) '5 


ORTH CPHin ,J) , Jr 1 ,191/0. , .355, .667, .95,1 .21 ,1 .471 ,1 .726,2,01 ,2,32 
,2.666,3.069,3.511,4.086,4.724,5,458,6.307,7.275,8,391,10.5/ 


DRTR CPHIC2,J) ,Jrl ,191/0. , .139, .291 , .451 , .63, .82,1 .031 ,1 .251 ,1 .479 
,1.721,1.97,2.224,2.478,2.77,3.055,3.391,3.72,4.12,4.49/ 


DRTfl tPHIC3,Jl ,Jrl ,19) /O. , .133, .271 , .42, .58, .759, .952,1 .165,1 .406, 
1.574,1.97.2.302,2.664,3.073,3.495,4,005.4,56.5.114,6,/ 


ORTH [PHIC4,J1 ,Jrl ,191/0. , .34, .51 , .659, .811 ,1 . ,1 .2,1 .41 ,1 ,639,1 .91 
1,2.224,2,556,2.903,3.297,3.760,4,26,4.77.5,3,6.2/ 


CONT. ON PG 
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DRTH CPHH5,J) ,Jrl .19)/Q, . .DIB, ,05B, .127, .215, .327, ,4BB, .649, .843, 
1 ,D7S.l ,344,1,644,1 .99,2.45,3,019.3.82,4,802,6,06,8,5/ 



CCNT. ON PG 3 

PG 2 OF 4 


JPL TECHNICAL REPORT 32-1523 









PB 3 DF 














PS5E I 


SUBROUfINE UBV(BN9PHA«RS«RG«V«BV9UBI 


SUBROUTINE Uavi3N«PHA«RS.R0>V»BW.UB) 


DIMENSION N3(<J2)»A« 3»M2I 

___A(J_iU_iAS_A5E__«Ii_iIL>_’_S 

DATA |AU»J)» J=l»<*2>/-26 

212,1 «-8.SS« 2.6Qi2.22« .72 
31.7C.2.6C. 1.3C»t.A9>3.6C 
AI2tJ)»S ARE 8-W»S 
__JDAJiA__(At2_» JJ_iJ.=J_iA2>_/ .6 3. 
l*87».83<.86>.80>.30f.8Q< 
.2_.7A»i_._2A»_.6A_»_.7i_»_.AO_. 
A13»J}»S ARE U-a*S 


«P-(I{6.19}»NP142} 


. 74 •-.42. -4, 41. -3.8 7», 2 l»-l, 7 3* 12, lei 3, 3s 
.-i_,.2.Ct6_.M.J_,C_C_.Sj,3CtlA,_l.JLi.._6j_lJ_eSjLlJ_»4_t_ 
•.aS..21.-l.ie«4.El»1.48*4.45>4.45«-?.lSe 
•-e.87.-1.16.4.CC.-l.Cl/ 

.»_eA7jL,62. ,2CKe22 »l^,3..3 

.80es80«.80«.80>.8Q.l.G4*.60t.62..73e.7l« 
_»„.A6 a. 6C . ,8j3 ._.£2_. .A5 a,AQa.AI jJl 


l.SOt.SSt. 40 ..4Qt.4C«.4Q< 

. _2_ s 7_5 j j_ 4_2 » ji_2 8_ » JL 4 D_» .AQjl > 2„8 . . 

NP*S ARE MODEL NUMBERS 
.__C AJ A__i NJPi_Ji, » J_= 1_» 4A1 /_C »_1 f 
«6>6.Et6.6«6.6f5.E<5.6iE« 


.4Q..40..4C«.40e.58. .30«.30e.34«.30e.35. 
.4J_....4C_._.2 5.,2 4j.,4A>_..2J.._.AC_.j,AC_.,^^^^ 


-2 »_3 »6 » 4 » 6.» 6.» 5 t.6»_l.»_Ej..6 J.6 »_6i 6_»6.i6_. 6_i j_6jf_ 
6f6t6t5i6.6.4/ 


DATA NB/0.- 100.-200.-30D.-301.-40C. -40 1.-402.- SCO. -501. -502.-503* 
Aji5il4 j_-5A§j.rACJ6trAL2 .-_5.6S_._-.6CA., -.5113 .1-5.11 ..-5.1 2_..rAQiLs-6C.lj_-6iL2>-f.Cijt_- 
*-604 1-605. -606.- 60 7. -60 8.-60 9. -61 0.-700. -701. -702. -703. -704. -70 5, 

.A-MCjLrACJL..:AC2 

PHI'S ARE -2. 5*L0GlC<PHi> ’S 

CAIA (P HI (i.Jl.J =1.19 1 /0. ..355. .66 7. .95. 1.21, 1.471.1. 726. 2 . 0 1. 2. 32 
«> 2.6 66. 3.06 9. 3.51 1 . 4.0 8 6. 4. 72 4.5.4 58.6. 30 7. 7.2 75. 8.391. 10.5/ 
.__E_A_l_A__iJ=_HXlA.U_ljA_=_ljJL9>./_0_,_._.133_.,A9X.j.45lj^O_t,,.82.._J_.E_lljLl*._25_lfJL._‘L2A 
*• 1.721. 1.97. 2. 224. 2. 478. 2. 7 7. 3.055.3.39 1. 3. 7 2. 4. If 2. 4.4 9/ 

__C_A_T A _J_F H_I l_3_._J_l_..U_=_l » I ? I/O . . .13 3 . , 21 1 . . 4_2 • _._5 6 . . 1 5 9 . ,_9 5_2 .J , 1 61»1. 4C_6jl 
* 1.674. i. 97. 2. ^0 2. 2,664. 3.073.3,495. 4.005. 4. 56. 5 . 1 1 4. 6./ 


*l .2, 224. 2,556. 2.903. 3.2 S7.3. 760.4.26.4, 77.5, 3.6.2/ 

_C_AJJ_lP_HI_lE_.Jl.J=1.16jyC.._.C 16r.C_56A.127..115._,J17_._.418_._,149_._^^^ 

*1,079. 1. 344. 1.644. 1 . 99 . 2 . 45 . 3 .0 1 9 • 3. 32 . 4 . 80 2 . 6,06. 8,5/ 

0 A1 A__i_P_H I_( A. J_) . jj =1 »_1 91/ 0 ... 2 8 . . 5A4 8_4 5_. l_._ 1 3 1 . I . _4_C 6 ._1_.6_S._2_._E.C_2_. 

• 2,34 8.2.742. 3. 1 1. 3. S2 2~. 3, 9^4 . 4, 52 4 . 5 . 1 1 4. 5.891. 7.160. 3.867. 12.0/ 

NBN = ( BN*1DC. 1-0.1 

DO 101 1=1.42 

If_{NBi_Il._N_E.MN_»__SO 1_0__ICC 

99 K = I 

S0_T_0__110 

100 IFa,EQ.4r} GO TO ^40 

101 CONTINUE . 

110 BV=A12.KI 

JUB_=_A_<_3_._K>_ 

NPFNrNPlKJ 

I_F(Mf_N._EQ_.C_J__6_Q__T0__2_l_3 

200 IPHA=PHA/10 
PHA=I 
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SUBROUTINE UB i BN « PH RS > R0« V « SV»U 81 


CIFF=TPHA-NPHa 

imELHi*.! 

WftL=DIFF*tPHI CNPFN tl*l l-PH I CNPF N . I > ) ■» PHI ( N FFN • I > 

&0^ASi.Z2£L 

211 V = An.U*S,>»ftLOElClRO> 

AO.AA.Zm 

22E W=A(l,KJ *5.*AL0G1CIRS*R0> *WAL 

S O TC 2 3C 

2MC PRINT 241 

— _2_4l__FOJ?fLA_lX2JiyillLLEJiAL__5LOD.T_JiUilBiRX 

B» = C. 

jia=iLc 

W = 99. 
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